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Abstract. I analyze the stochastic effects introduced by the sampling of the stellar initial mass 
function (SIMF) in the derivation of the individual masses and the cluster mass function (CMF) 
^jry from broad-band visible-NIR unresolved photometry. The classical method of using unweighted 
UBV photometry to simultaneously establish ages and extinctions of stellar clusters is found 
, to be unreliable for clusters older than w 30 Ma, even for relatively large cluster masses. On 
the other hand, augmenting the filter set to include longer-wavelength filters and using weights 
for each filter increases the range of masses and ages that can be accurately measured with 
unresolved photometry. Nevertheless, a relatively large range of masses and ages is found to be 
dominated by SIMF sampling effects that render the observed masses useless, even when using 
UBVRIJHK photometry. 
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1. Description 

This work is the third of a series in which we are analyzing the possible biases present 
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in mass functions (Mai'z Apellaniz & Ubeda 2005 |Maiz Apellaniz 2008[ ). I have used a 



combination of analytical approximations and Monte Carlo simulations to study the effect 
of the stochastic sampling of the stellar IMF (SIMF, assumed to be Kroupa) in the deter- 

• mination of masses and ages of unresolved stellar custers from broad-band photometry. 
For this purpose I have used the new (3.1) version of CHORIZOS ( |Maiz Apellaniz 2004P , 
which incorporates an evolutionary synthesis module. For a given mass and age, I gen- 
erated a minimum of 10 000 realizations of the SIMF using solar metallicity [a] Geneva 

• isochrones for massive stars and [b] Padova isochrones with new AGB treatment for low 
and intermediate masses. In each case I measured the mass (and age, if appropriate) 
assuming that the cluster has a well-sampled SIMF. From the individual realizations I 
derived [a] the observed cluster mass and age distributions for clusters of a (fixed) real 
mass and age and [b] the observed cluster mass functions (CMFs) for a real truncated 
power-law CMF with a slope of 7 = —2.0. Three different cases in order of complexity 
were considered: [a] single-filter observations of clusters with known age and extinction, 
[b] multi-filter observations of clusters with unknown age and known extinction, and [c] 
multi-filter observations of clusters with unknown age and extinction. In the first case 
masses were computed directly by converting filter-convolved luminosities to masses while 
in the second and third cases a \ 2 minimization code (CHORIZOS) was used to derive 



the masses and ages. The reader is referred to Mai'z Apellaniz (2009) for further details 
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Figure 1. Results of the Monte Carlo simulations of fM \M(mo\m) and Gaussian and Gamma 
fits for 10 Ma clusters observed in V with four real masses: 300 M© (upper left), 1000 Mq (upper 
right), 10 000 M (lower left), and 100 000 M (lower right). The Poisson mass 9 for this filter 
and age is 566 Mq. For large real masses, both fits provide reasonable approximations but for 
small values of m, the Gamma approximation is significantly better. The vertical dotted lines 
show the values of the real mass m and, when visible, of 8 and tti ,lll (the observed cluster 
mass obtained when one mistakes the brightest possible star in the isochrone for a cluster). 

2. Case 1: Single filter, known age and extinction 

• For a given cluster mass (m), the observed cluster mass (m Q ) distribution can be 
approximated by a Gamma function better than by a Gaussian distribution (Figure [I}: 

fM \M(m \m) = Am^e-^e-^ 6 . (2.1) 

• The Poisson mass 9 is a measurement of the stochasticity and (strongly) depends 
on the age and the filter used to derive m D (Table [2]) . 

m — m ; <r mo = V 9m. (2-2) 

• The large differences in 9 are a consequence of the different fraction of stars located 
above the median mass for the luminosity in each filter for a given isochrone (Figure (2) ■ 

• The observed CMF, /m ( to o)i differs from the real (truncated power-law) CMF, 



S266. The validity of observed cluster masses and ages 



3 



3-10° 



1-10" =■ 
3-10 5 — 
MO 5 — 
3-10 4 — 
MO 4 =■ 
3-10 3 — 
MO 3 — 
3-10 2 — 
MO 2 — 
3-10' — 
MO' =- 
3-10° — 
1-10° =■ 
3-1 0" 1 — 



10" a 
10 80 a 




. 10 lou a-AGB6 

fa U median 
-fa B median 
V median 
-fa R median 
/ median 
fa /median 
fa //median 
fa ^median 

1 in 500 

1 in 1000 

50000 



20000 



10000 



5000 



T e „ (K) 



Figure 2. Isochrones for seven of the nine ages used in this work. For each isochrone, symbols 
of different colors are used to indicate the median mass for the luminosity in the eight filters 
U, B, V, R, I, J, H, and K assuming a Kroupa IMF between 0.1 M Q and 120 M . The two 
dashed lines join the points in each isochrone beyond which 0.2% and 0.1%, respectively, of the 
remaining stars in a well-sampled Kroupa IMF are located. The two remaining isochrones are 
not included for the sake of clarity. 



/iW-(m), in two aspects (Figure[3]): [a] For large masses, /m (w ) is similar to /m(to) but, 
as the cluster mass decreases, at one point a "hump" or overdensity appears, [b] To the 
left of the hump, /m„ {m ) asymptotically approaches a power law with a slope larger 
than —1.0, with some clusters having m values lower than the real cutoff in to. 

• Using an analytical approximation, it can be shown show that if one is willing to 
tolerate a systematic error in the CMF slope 7 up to A7, then one can only analyze 
observed masses: 

mo > 0.57(7+ 1)0/ A7, (2.3) 
e.g. for 7 = -2.0 and A7 = 0.1, to q > 10 8. 
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Table 1. Values of 6 (in M©) for different ages and niters for clusters observed with a single 
filter and known age and extinction. The labels below each age indicate the type of the brightest 
possible star. For a given filter, 9 starts at a low value in the MS phase, rapidly grows in the 
LBV phase, decreases during the RSG and early AGB phases, and then experiences a moderate 
growth when the cluster becomes old. For a given age, 9 increases as a function of wavelength 
with the only exception of the MS phase, where the effect is the opposite (but rather weak). 
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Figure 3. Observed (continuous lines) and real (dotted lines) CMFs for the minimum and 
maximum values of 9 in Table [2] The left panel corresponds to 1 Ga clusters observed with U 
and the right panel to 316 Ma clusters observed with K. Each panel shows four cases of the 
lower mass cutoff for /Af(m) from 1 Mq to 1000 M©. As 9 increases, the hump becomes more 
pronounced and moves towards the right. 

3. Case 2: Multiple filters, unknown age, and known extinction 

I used three execution types: 

1. UBVRIJHK photometry with single-band ^-dependent weights for each filter. 

2. UBVRIJHK photometry with constant weights. 

3. UBV photometry with single-band ^-dependent weights for each filter. 

The observed age distributions for 10 4 M clusters are shown in Figure 2] The main 
results are: 

• The ages calculated with UBVRIJHK photometry and constant weights (i.e. the 
standard method) are highly uncertain even for 10 4 Mq clusters. 

• Using weights based on the stochasticity of each filter significantly narrows the 
observed-age distribution for 10 4 M Q clusters, thus allowing for the obtention of ages 
with relatively small uncertainties. 
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Figure 4. Distribution of observed ages for the nine input ages and three execution types (1: 
upper left, 2: upper right, 3: bottom) for clusters of 10 4 Mq and unknown ages. At the top of 
each panel the dashed line indicates the real age while the symbols and error bars provide the 
median and inferior and superior uncertainties (1-sigma equivalents) of each distribution. 



• UBV and UBVRIJHK photometry provide relatively similar results, with the 
former having the advantage for old clusters and the latter for intermediate-age ones. 

• The AGB1 stage (31.6 Ma) has significantly larger uncertainties than the rest. 



4. Case 3: Multiple niters, unknown age and extinction 

I used two execution types: 

1. UBVRIJHK photometry with single-band ^-dependent weights for each filter. 

2. UBV photometry with single-band ^-dependent weights for each filter. 

The observed age distributions for 10 4 M Q and 10 5 M Q clusters are shown in Figure [51 
The main results are: 

• As opposed to the cases where extinction is known, the addition of RIJHK pho- 
tometry provides a significant improvement in the observed ages. 

• U BV photometry alone is not sufficient to accurately determine the ages of unre- 
solved stellar clusters, even for massive ones, if extinction is unknown. 

• For some ages (1-10 Ma, 100 Ma, 10 Ga), UBVRIJHK photometry provides rela- 
tively accurate observed ages for m = 10 4 M Q . For the rest, higher masses are needed. 
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Figure 5. Distribution of observed ages for the nine input ages and two execution types (1: left 
plots, 2: right plots) for clusters of 10 4 Mq (top) and 10 5 Mq (bottom) and unknown ages and 
extinctions. At the top of each panel the dashed line indicates the real age while the symbols 
and error bars provide the median and inferior and superior uncertainties (1-sigma equivalents) 
of each distribution. 

5. Conclusions 

• SIMF sampling effects can introduce large biases in the determination of masses and 
ages from unresolved photometry. 

• One should use weights based on the relative stochasticity of each filter when using 
X 2 -minimization techniques to calculate cluster properties from unresolved photometry. 

• The addition of a broad-photometric baseline (U to K) significantly reduces uncer- 
tainties when extinction is determined from the data. 

• For clusters younger than 100 Ma, a critical (post-RSG) stage exists around 30 Ma 
where ages are especially difficult to determine. 
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